Absorbing boundaries in numerical solutions of the time-dependent Schrodinger 
equation on a grid using exterior complex scaling 
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We study the suppression of reflections in the numerical simulation of the time-dependent 
Schrodinger equation for strong-field problems on a grid using exterior complex scaling (ECS) as 
absorbing boundary condition. It is shown that the ECS method can be applied in both the length 
and the velocity gauge as long as appropriate restrictions are applied in the ECS transformation of 
the electron-field coupling. It is found that the ECS method improves the suppression of reflection 
as compared to the conventional masking technique in typical simulations of atoms exposed to an 
intense laser pulse. Finally, we demonstrate the advantage of the ECS technique to avoid unphysical 
artifacts in the evaluation of high harmonic spectra. 

PACS numbers: 32.80.Rm, 33.80.Rv, 42.50.Hz 
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I. INTRODUCTION 

With the progress in laser technology in recent years 
P, [1] , the focused laser field strengths increased rapidly 
to exceed the strength of the Coulomb fields that bind 
the electrons in the ground state of an atom or molecule. 
Such intense light fields have led to the discovery of 
novel aspects of light-matter interactions, such as multi- 
photon ionization, above-threshold ionization, high har- 
monic generation, multiple ionization to high charge 
states, Coulomb explosion etc. The lowest-order per- 
turbation theory of light-matter interaction is known to 
break down at light intensities above about / ~ 5 x 10^^ 
W/cm^ for near optical wavelengths 0]. Several non- 
perturbative theories of laser-matter interaction, such as 
the numerical solution of the time-dependent Schrodinger 
equation (TDSE) 0, i, %M, ab-initio methods based 
on the Floquet theorem [1, [3, basic set expansion 
methods [ i j] . [l^. [Tsj or iS-matrix and related theories 
[3 [S Il6l Il7l|. have been developed. Among these the 
solution of the TDSE on a time-space grid is consid- 
ered as a rigorous and powerful approach. For the in- 
vestigation of most of the above-mentioned intense-field 
phenomena the simulation of an atom with N electrons, 
which would require the solution of a set of 3 x iV di- 
mensional partial differential equations, can be well ap- 
proximated using the singlc-active-electron (SAE) model 
[J|. In the latter only one of the electrons of the atom 
is considered to become active and to interact with the 
external field, which reduces the numerical problem to 
at most three dimensions. Any symmetry of the external 
field may further reduce the dimensionality. 

Accurate solutions of time-dependent problems on the 
grid require not only relatively dense grid points but, 
in general, also a huge spatial extension of the grid to 
account for the release and the motion of the electron 
sub-wave packets in the field. These factors result in 
large memory and CPU requirements for the numerical 
solution at high field intensities. Fortunately, the im- 
portant dynamics of several intense-field processes occur 
on a spatial volume close to the atom or molecule. A 



few examples are the determination of ionization rates 
[1] and high harmonic spectra P| or the identification of 
dominant pathways to single [18] or nonsequential dou- 
ble ionization 19] of a molecule. These studies can be 
therefore performed on a relatively small grid neglecting 
the exact form of the outgoing ionizing parts of the wave 
function. It however requires to suppress refiections from 
the edges of the numerical grid, which can cause artifi- 
cial effects, e.g. in the form of spurious harmonics in high 
harmonic spectra 

In calculations based on the numerical propagation of 
wave packets in different areas of physics and chemistry 
several techniques have been proposed to compensate 
for reflections, including masking functions (or equiva- 
lently absorbing imaginary potentials) d, [l^, [2l| , repeti- 
tive projection and Siegert state expansion methods [22[, 
com plex coordinate rotation or exterior complex scal- 
ing |23l . [13] and others. In the numerical solution of 
the time-dependent Schrodinger equation of atoms in 
intense laser fields masking functions or absorbing po- 
tentials are the most commonly used techniques (e.g. 
US H, 113, m, HI). The ECS technique is rarely 
used in this context up to now [H, HI, Q , which might 
be due to the fact, that its application in the electric field 
(or length) gauge has been thought to be not fruitful [s^ . 

In this paper we re-examine the implementation of the 
exterior complex scaling (ECS) method [13, [Ml as ab- 
sorbing boundary condition in simulations of strong-field 
problems. In particular, we focus on the application of 
the ECS technique in different gauges, namely the length 
and the velocity gauge. On the basis of results of simula- 
tions of ID model atom exposed to an oscillating external 
field, we investigate appropriate restrictions of the ECS 
transformation of the external-field coupling to avoid un- 
desired effects in the absorbing area. Results for the 
probability density and momentum distributions for the 
interaction of the hydrogen atom with an intense laser 
field are then compared with those obtained using the 
standard masking function technique. Finally, we con- 
sider a typical strong-field problem by using the ECS 
method to calculate high harmonic spectra. 
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FIG. 1: Scheme of the ECS coordinate transformation. Real 
coordinates are used in the interior box, defined by zi < z < 
zi and p < pa- The zones in gray mark the areas, where one 
or both coordinates are complex. 



II. EXTERIOR COMPLEX SCALING (ECS) AS 
ABSORBING BOUNDARY 



The complex scaling method has been widely used in 
physics and chemistry (for a review, see (39|). e.g. in 
the theory of potential scattering [4ffl , calculation of res- 
onances in atoms and molecules |4l| or the calculation 
of cross sections in scattering processes [i^. According 
to this method the radial coordinate of the particles are 
scaled by a complex phase factor, which distorts the spec- 
trum of the Hamiltonian such that the continuous spec- 
trum is rotated in the complex energy plane and the dis- 
crete resonance eigenvalues are revealed. For our aim to 
introduce absorbing boundaries at the edges of a numeri- 
cal grid in time-dependent simulations we make use of an 
extension of the complex scaling method, iiamely the ex- 
terior complex scaling (ECS) technique [31|, in which the 
spatial coordinates are only scaled outside some distance 
from the origin. 

As discussed at the outset the ECS technique has been 
used before [13, [s^l in time-dependent studies of 
electron impact ionization as well as of the motion of 
a charged particle in dc and ac electromagnetic fields. 
In the latter context it has been mentioned [s^ that it 
appears to be not fruitful to apply ECS in the electric 
field (or length) gauge. We now revisit this question and 
analyze below how ECS can be used in both the length 



and the velocity gauges in simulations of atoms exposed 
to intense laser pulses. 



A. Implementation 

First, we illustrate the implementation of the ECS 
technique in ID and 2D time-dependent calculations, it 
is straightforward to extend it to higher dimensions (c.f. 
[U, [3^ . Let us consider the non-relativistic electron dy- 
namics in a time-independent (Coulomb) potential Vq 
and an external field governed by the time-dependent 
Schrodinger equation in cylinder coordinates as (Hartree 
atomic units, m,, = h = e = I, are used throughout): 

. d 



Vitmz,p;t) (1) 
with the time-indepentent Hamiltonian 

Ho{z,p) 



1 Id 1 92 , 



2dz^ 

and time-dependent external-field coupling, given in 
length or velocity gauge, as 



length gauge 
velocity gauge 



(3) 



Here E{t) is the electric field and A{t) is the vector po- 
tential of the external electromagnetic pulse linearly po- 
larized along the z-direction. 

The ECS transformation on the two coordinates z and 
p can be given by (c.f. Fig. [J): 



zi + {z — zi) exp(i?7) as 2; < zi 

z as zi < z < Z2 

Z2 + (z - Z2) exp(i77), as z > Z2 

p as p< pq 

Po + (p - Po)exp(iry), as p > po 



(4) 
(5) 



where 77 is the scaling angle with < 77 < 7r/2. zi, Z2 and 
po are labeled in Fig. [ijand define the size of the interior 
box (zi < z < Z2 and p < po) within which both spatial 
coordinates are real. Outside (gray zones in Fig. [1]) one 
or both coordinates are complex. 

It is the aim of the ECS method to transform the out- 
going wave into a function, which falls off exponentially 
outside the interior box, while the wave-function keeps 
unchanged in the region where the coordinates are real 
^A, 36, J57.] . In case of the present problem of an atom 
exposed to an oscillating linearly polarized electromag- 
netic pulse, we therefore investigate whether or not the 
transformed solution of Eq. [1] shows this desired behav- 
ior. To this end and without loss of generality, we restrict 
our analysis to the Z direction, i.e. the direction of the 
external field. The time-dependent solution of the ID 
analogous of Eq. [1] in the complex area can be written 
as: 

ij{t + At) - exp[-i{Ho{Z) + V(t))At]i;{t). (6) 



The time-independent operator in Eq. ([6]) is given by: 



exp[-ii7o(^)Ai] 
exp 



.cos(2?7) 
« 7i Pz 



exp 



sin(277) 



exp[-iRe{Vo)At] exp[Im{Vo)At] 



(7) 



As discussed by McCurdy et al. 33] the exponent in the 
second factor on the right hand side of Eq. ([7]) is always 
negative if < 77 < 7r/2 and provides aheady the desired 
decay term. It is therefore important to note that the 
wave- function wiU be basically absorbed in the complex 
area due to the transformed kinetic operator term as long 
as there are no counteracting effects from other terms in 
the Hamiltonian. In general, it is therefore required that 
Jm(Vb) < such that the last term in Eq. (O acts 
as an absorbing potential. In practice, the discontinuity 
in the real and imaginary part of the potential terms 
introduced by the complex scaling factor can generate 
some small reflections. An efhcient way to avoid this 
numerical problem appears to keep the Z-coordinate in 
the potential term untransformed as a real number in the 
absorption area. Please note that it is unproblematic to 
abandon the additional absorption effect of the imaginary 
potential, since the decay of the wavefunction is already 
ensured via the kinetic operator term. 

The time-dependent potential in Eq. ([6]) can be writ- 
ten as 

exp[~iE{t)ZAt] 

= exp[~iE{t){zi^2 + {z - -21,2) exp(i77))Ai] 
= exp[E{t){z — zi_2) sin?7At] 

X exp[-~iE{t){zi^2 + {z ~ 21^2) COST]) At] (8) 



in the length gauge or as 
/A{t) 



exp 



V c 



exp(-i?7)p^ At 



exp 



sm?7 PzAt 



exp 



-I cos 77 PzAt 

c 



(9) 



in the velocity gauge. In the right sides of Eqs. ^ and 
(ini), the second factors are oscillatory ones, and are sim- 
ple scaled versions of the external-field coupling as the 
wavefunction enters the complex area. The first factors, 
however, can act both as an absorber or as an undesired 
source, depending on the sign of the exponent. It is de- 
termined in the length gauge by the instantaneous sign of 
the oscillating electric field, while in the velocity gauge, 
it equals the sign of the product A{t)pz- Next, we will 
analyze the effects of the oscillating field coupling term 
in numerical simulations and show that an undesired ex- 
plosion of the wavefunction can be avoided by using the 
standard untransformed coupling even in the absorption 
region. 
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FIG. 2: (color online) Temporal evolution of the probability 
density from the numerical simulation of a ID model atom in- 
teracting with a high-frequency field (left hand column) and 
a Ti:sapphire field (right hand column). Shown is a compari- 
son of the results obtained using the external field coupling in 
the velocity gauge (upper row) and in the length gauge (lower 
row). 



B. ECS technique in the length and the velocity 
gauges 

In order to investigate the effects of the oscillating field 
on the transformed wavefunction we have performed ID 
test calculations with a soft atomic model potential 



1 



(10) 



using two different external fields, namely a high- 
frequency electric field given by 



^ _ ( El sm{uJit)t/5T, as t<5T 
^ \ El sincjii, as t > ST 



(11) 



where Ei ~ 0.5, cji = 0.5 and T — 27r/a;i, and a 3-cycle 
low- frequency Ti: sapphire laser pulse 



E ^ E2 COs{uJ2t) SIT? {Tit /L2) 



(12) 



with E2 = 0.1 a.u., LU2 — 0.057 a.u. and L — 330 a.u.. 
The respective vector potentials are derived from the 
electric field expressions in Eqs. pT|) and ([12]). Please 
note that the former field, which is smoothly turned on 
over 5 optical cycles, is similar to the field form used by 
McCurdy et al. 33]. The real part of the calculation box 
is restricted in both cases by zi = —Z2 — —25 a.u., with 
the complex part extending over 12.5 a.u. on both sides 
of the grid. 

In Fig. [2] we present the temporal evolution of the 
electron density distributions in the high-frequency field 
(left hand column) and in the Ti:sapphire field (right 
hand column). The panels in the upper and the lower 
row show the numerical results obtained in the length 



4 




100 200 300 

t(a.u.) 



1D0 2D0 300 

t(a.u.) 





100 200 

t(a.u.) 



100 200 

t(a.u.) 



FIG. 3: (color online) Same as Fig. [2] but using the untrans- 
formed standard field coupling on the entire grid. 



and the velocity gauge, respectively. The effect of the 
first factor in Eqs. ([8]) and ([9]) is most clearly seen in 
the results of the low-frequency calculations, where we 
observe an explosion of the wave-funtion in the second 
part of the evolution after a significant part of the wave- 
function has entered the complex area. This is obviously 
due to the fact that the term acts as a source over a half 
cycle of the pulse. As expected above the same unphysi- 
cal feature is found in the length as well as in the velocity 
gauge. In the high-frequency case (left-hand column) the 
results obtained in the two gauges are again similar, but 
we do not observe a significant amount of explosion. This 
difference as compared to the Ti: sapphire calculations is 
probably due to two factors: First, in the high-frequency 
simulation the probability density, which enters the com- 
plex area, is smaller than in the low-frequency calcula- 
tions. Second, the rapid change of the sign of the electric 
field or the vector potential may effectively prevent an 
explosion, since the complex factor quickly changes be- 
tween decay and source nature. 

From the results presented above we may therefore in- 
fer that the ECS transformation may lead to unphysical 
results in both gauges due to the oscillating nature of the 
external field. The most straightforward strategy to cir- 
cumvent this problem appears to simply neglect the ECS 
transformation of the coordinates in the field coupling 
term. This means in other words that the untransformed 
standard field coupling (c.f. Eq. ([3])) is used over the en- 
tire grid including the absorbing area. Using this strategy 
there is no risk to create a source term while the desired 
decay of the wavefunction should be still achieved via the 
kinetic operator term. 

In order to test our expectations we have repeated the 
simulations by dropping the term exp{irj) in Eq. ([H]) and 
exp(— 177) in Eq. ([9]), the results are shown in Fig. [S] The 
comparison with the respective panels of Fig. [3] show im- 
mediately that the desired effect is achieved. In partic- 
ular in the low-frequency case (right hand column) the 



wavefunction is absorbed at the edges of the grid and no 
signature of explosion is seen anymore. 

We may therefore conclude that the exterior complex 
scaling technique can be applied as an absorbing bound- 
ary to time-dependent simulations on laser-atom simu- 
lations in the length as well as in the velocity gauge as 
long as the complex factors in Eqs. (O and Q are re- 
moved. We may note that McCurdy et al. [33| reached 
to a similar conclusion for the velocity gauge as they did 
not transform the momentum pz to the complex plane. 
Our analysis above shows that an analogous restriction 
is possible in the length gauge as well. Our test calcula- 
tions have shown that, in general, after omission of the 
unstable factors calculations in the length gauge show 
slightly better results than those in the velocity gauge. 
We therefore restrict ourselves below to the length gauge 
only. 

Thus, the above implementation of the ECS technique 
coincides with the desired absorbing boundary condition 
in time-dependent strong-field calculations. It has been 
shown in the application of ECS to the time-independent 
Schrodinger equation (e.g. [H, and for review [24]) 
that using a sharp exterior scaling the derivative discon- 
tinuity at the boundary is handled exactly, as long as 
the boundary is chosen to coincide with a grid point. We 
have adapted this strategy in the time-dependent calcula- 
tions. In test calculations we have found that the absorp- 
tion effect at the boundaries is almost independent of the 
scaling angle 77, in the present calculations we have used 
rj = jZ. Finally, before proceeding with a comparison 
of the ECS results with those obtained using the stan- 
dard masking function technique we may note that spe- 
cial care has to be taken in the representation of the first 
and second derivatives in the transformed Schrodinger 
equation (for a detailed discussion, see [111), which we 
have approximated using Lagrange interpolating polyno- 
mials 14311. 



C. Comparison of ECS and masking function 
techniques as absorber 

In order to analyze the efficiency of the ECS method as 
absorber in time-dependent simulations as compared to 
the masking function technique we present in this section 
results of calculations for the interaction of the hydrogen 
atom with a linearly polarized laser field. We compare 
the results obtained using the ECS method with those ob- 
tained using conventional masking functions of the form: 



M 



cos* 



\X - Xi\lT 



(13) 



where d is the length of the absorbing region, over which 
M changes smoothly from 1 to 0, and Xi is the boundary 
point. Such a function has been applied at all boundaries 
of the numerical grid (c.f. gray zones in Fig. [T]). In 
the course of the calculations we have tested different 
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FIG. 4: (color online) Logarithmic contour plots of the elec- 
tron density distributions from the numerical solution of the 
2D TDSE of the hydrogen atom exposed to an intense laser 
pulse ai t = 165 a.u. (upper row) and t = 330 a.u. (lower 
row) . Results obtained using ECS (left hand panels) are com- 
pared with those obtained using the masking function (right 
hand panels) as absorber. Field parameters as in Eq. (I12|l . 



masking functions, the results presented below are found 
to be rather insensitive on the form of masking functions. 

We use the time-dependent Schrodinger equation in 
the length gauge and the Coulomb potential of the hy- 
drogen atom, 

V{z,p) = E{t)z- J-^ . (14) 

The electric field profile is given by Eq. p^ . The field 
parameters are the same as before. The grid parameters 
are Az — 0.1 a.u., Ap = 0.2 a.u. and the time step At = 
0.1 a.u.. The initial ground states have been obtained 
via imaginary time propagation [i^ . The absorber is 
applied at po = 22 and zi = —Z2 = —10. The width of 
the absorber is chosen to be 20% of the grid size. 

To demonstrate that the reduction of reflections is 
present in the solution of the TDSE, we present in Fig. 
|4] probability density distributions. From the compari- 
son between the results for ECS (left column) and for 
the masking function (right column) at t = 165 a.u. and 
t = 330 a.u. in the upper and lower row, respectively, 
the difference in the efficiency of the two absorbers is 
clearly visible. The distributions obtained with the mask- 
ing function show interference patterns due to the reflec- 
tions at the boundaries, which are not seen in the results 
for the ECS absorber. Note that at i = 165 a.u. the 
wave packet has reached (and is reflected from) the up- 
per boundary in ^-direction only, while at the end of the 
pulse reflections in all directions have appeared. 
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FIG. 5: (color online) High harmonic spectra from the inter- 
action of the hydrogen atom with an intense linearly polar- 
ized 3-cycle laser pulse at Eo = 0.1 a.u. and cu = 0.057 a.u.. 
Shown is a comparison between the results obtained from the 
(unrestricted) reference calculation (black solid line) and cal- 
culations using the ECS (red dashed line) and the masking 
function (blue dotted line) as absorbers at zi = —Z2 = —80 
a.u.. In transversal direction the grid was large enough to 
avoid additional reflections from this boundary, (a) Full spec- 
trum and (b) cut-off region. 



III. CALCULATION OF HIGH HARMONIC 
SPECTRA 

Finally, we apply the ECS absorber to a typical 
intense-field phenomenon, namely the evaluation of high 
harmonic spectra. High harmonic generation (HHG) 
is an important process for laser frequency conversion 
and the generation of attosecond pulses (for reviews, see 
e.g. (isl. |46||). According to the semiclassical three-step 
rescattering picture (47l . [48l . and confirmed by quantum- 
mechanical calculations [4!l| , HHG can be understood as 
the ionization of an electron by tunneling through the 
barrier of the combined Coulomb and laser fields, fol- 
lowed by the acceleration of the electron in the field, 
which may cause, for linear polarization of the field, a 
return of the electron to and its recombination with the 
parent ion under the emission of a harmonic photon. On 
the basis of this picture it is reasonable to limit the grid 
size of an ab-initio calculation of high harmonic spectra 
via the time-dependent Schrodinger equation, since be- 
yond a certain distance from the nucleus outgoing wave 
packets are expected to have no effect on the high har- 
monic spectra. 

We have performed simulations of the 2D TDSE for 
the hydrogen atom in an intense linearly polarized laser 
pulse, given by the pulse form in Eq. with Eq — 0.1 
a.u. and lu — 0.057 a.u., and evaluated high harmonic 
spectra as the Fourier transform of the time-dependent 
dipole moment. Note that the dipole moment has been 
determined over the interior box (i.e. without the ab- 
sorbing regions) . In order to analyze the effect of refiec- 
tions from the edges of the grid along the polarization 
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direction on the spectra, we compare in Fig. [5] the re- 
sults of three simulations, namely the full calculation as 
a reference (black solid line), performed on a sufficiently 
large simulation box to prevent reflections at the bound- 
aries, and calculations using ECS (red dashed line) and 
masking function (blue dotted line) as absorbers. The 
calculations with absorbers have been restricted along 
the polarization axis by choosing zi = —zi = —80 a.u., 
which exceeds the maximum excursion of the classical 
electron trajectories of 63.6 a.u.. The absorbing part of 
the grid has been chosen to extend over an additional 20 
a.u. at both ends. In the transversal direction the grid 
was chosen large enough to avoid reflections from this 
boundary. 

The results in Fig. [5^) show the typical high har- 
monic spectrum with a plateau and a cut-off at A'^ = 
{Ip + 3.17J7p)/w « 51, where Ip — 0.5 is the ionization 
potential of the hydrogen atom and Up = Iq/Au)"^ = 0.77 
is the ponderomotive potential. From the comparison in 
Fig. [5] it is seen that the harmonics in the plateau do 
not differ significantly. There are small deviations in the 
minima between the harmonics obtained from the simu- 
lation with the masking function but the maxima appear 
to be unchanged. The effects of the reflections become 
visible at and beyond the cut-off, this region is enlarged 
in Fig. ^p). While the results from the ECS calculation 
almost agree with those from the full calculation over a 
decrease in the signal of two orders of magnitude, the re- 
sults evaluated with the masking function start to differ 
near the cut-off and the deviations increase up to an order 
of magnitude in the signal beyond the cut-off. This "ar- 
tificial" increase of the HHG signal in the cut-off region 
results from those parts of the wavcfunction reflected at 
the boundary, which return to the nucleus and give rise 
to harmonics without physical meaning. Note that an ac- 
curate calculation of the harmonics in the cut-off regime 
is e.g. important for analysis of the generation of single 
attosecond pulses (e.g. |50|). 

The effect of the reflections even increases for longer 
pulses, in which several wave packets reach the bound- 
aries. This is seen from the results, presented in Fig. [SI 
where harmonic spectra obtained for laser flelds with a 
constant envelope, E = Eq sin(ci;t) with Eq = 0.05 and 
Lu — 0.057, having 3 (upper row), 6 (middle row) and 
12 (bottom row) cycles are shown. Results obtained us- 
ing ECS and masking function as absorber techniques 
are shown in the left and right hand panels, respectively. 
The absorbers were placed at zi = —Z2 = —35 a.u. (here 
the maximum excursion length is 30.8 a.u.), with an ad- 
ditional absorbing region of 6 a.u. The grid size in ^- 
direction was big enough not to influence the results. 
While the results obtained for the 3-cycle fleld show the 
expected cut-off at the 19th harmonic, spurious harmon- 
ics due to reflections appear beyond the cut-off and in- 
crease in signal as the number of fleld cycles increases. 
The comparison shows that for the ECS results the con- 
trast ratio between the false and the plateau harmonics is 
about lO"'^ in all cases, while for the masking functions 
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FIG. 6: Comparison of harmonic spectra calculated using the 
ECS (left column) and the masking function (right column) 
teclmique in laser pulse with constant amplitude having 3 
(upper row), 6 (middle row) and 12 (lower row) cycles. 



it increases, especially for the highest false harmonics, 
giving the impression of a second unphysical plateau. 



IV. SUMMARY 

In summary, we have investigated the implementa- 
tion of the exterior complex scaling technique as an ab- 
sorber in the numerical solution of the time-dependent 
Schrodinger equation for strong-fleld problems on a grid. 
Our analysis has shown that the ECS technique can be 
applied in both the length and the velocity gauge as long 
as the untransformed field coupling is used on the entire 
grid including the absorbing area. It is found that the 
decay due to the ECS transformation in the kinetic op- 
erator term is sufficient to efficiently reduce reflections 
at the grid boundaries. A comparative study has shown 
that in this implementation a significantly better sup- 
pression of refiections can be achieved as using the con- 
ventional masking function method. By application of 
the ECS method to the evaluation of high harmonic spec- 
tra differences in the suppression of artifacts, e.g. in form 
of spurious harmonics, is demonstrated. The simple test 
cases considered here should capture the essence of the re- 
fiection problem, exterior complex scaling in both length 
and velocity gauge can be therefore considered as an effi- 
cient absorption technique for numerical time-dependent 
solutions in higher dimensions and/or of more complex 
processes too. 
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